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Abstract 

Here, it is found that dark energy and dark matter emerge from 
the gravitational sector, if non-linear term of scalar curvature is added 
to Einstein-Hilbert lagrangian. An equation of state for dark energy, 
having the form p de = — p de + f(a) (with Pde(Pde) being the pressure 
(density) for dark energy, f(a) being a function of scale factor a(t) and 
t being the cosmic time) is explored. Interestingly, this equation of 
state leads to a phantom barrier wa e = Pde/ Pde = — 1 at a = a w . It is 
found that when a < a w ,Wd e > — 1 and w dc < — 1 for a > a w , showing 
a transition from non-phantom to phantom phase at a = a w < a^a^ 
being the current scale factor of the universe). 
PACS no. : 98.80.Cq 
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Late cosmic accleration is the most remarkable astronomical obser- 
vation in the recent past [1, 2]. Theoretically, it is possible when the 
universe is driven by a fluid having negative pressure and equation of 
state parameter (EOSP) w dc < —1/3. Observational data also indicate 
that Wde should be close to —1 (the phantom barrier) and most proba- 
bly less than — 1 for the current universe. The source of this mysterious 
fluid is still in dark. So, many phenomenological models, suggesting 
possible sources of dark energy (DE), have been proposed in the past 
few years [3, 4]. In spite of non-gravitational sources of DE, non-linear 
terms of curvature are also proposed as gravitational alternative for 
DE [5]. A recent comprehensive review on DE is available in [6]. 

In [5] , non-linear terms of curvature are taken as DE lagrangian a pri- 
ori and its various consequences are discussed. In the present model, the 
story is different, Here also, investigations begin from the modified 
gravity stemming from addition of non-linear term to Einstein-Hilbert 
lagrangian, but unlike [5], here, non-linear term of curvature is not 
taken as a DE source. Rather, DE emerges spontaneously as a com- 
bined effect of linear as well as non-linear terms of Ricci scalar cur- 
vature [7, 8], whereas, in [5], only non-linear term of curvature con- 
tributes to DE. Interestingly, here, it is found that curvature can be a 
possible source of dark matter too. In [9], various equations of state 
(EOS) for DE, dependent on Hubble's expansion rate H and its deriva- 
tives, are proposed. Contrary to this, in what follows, EOS for DE is 
not proposed but derived. EOS for DE, obtained here, has the form 
Pdc = -Pdc + f{a) ( with p de ,Pde, a(t) and t being dark energy density, 
pressure, scale factor and cosmic time respectively). It is interesting to 
see that a = a w < ao (ao is the current scale factor) gives the phantom 
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barrier Wd c = — 1. It is found that when a < a w , — 1 < Wd c < — 1/3 
and Wdc < — 1 for a > a w . It shows a transition from a quintessence 
to phantom phase at a = a w . In [10], a possibility of this type of tran- 
sition is discussed taking EOS in Jordan frame. The present work is 
different from [10] as no EOS is taken here a priori, rather it is de- 
rived yielding phantom barrier and quintessence to phantom transition 
spontaneously. 

Natural units (h — c — 1) are used here with GeV as the fundamental 
unit, where h and c have their usual meaning. In this unit, it is found 
that lGeV" 1 = 6.58 x 10- 25 sec. 

2. The action for higher-derivative gravity is taken as 

S 9 = J d 4 xV^[-J^ ~ aR(2+r) l W 
where R is the Ricci scalar curvature, G = Mp 2 (M P = 10 19 GeV is 
the Planck mass). Moreover, a is a coupling constant having dimen- 
sion (mass)~ 2r with r being a positive real number. As (2 + r) > 0, 
instability problem does not arise like the model containing the R^ 1 
term [11]. But, as mentioned above, here approach is different from 
the papers [10, 11]. 

The action (1) yields gravitational field equations 

-^ (2+r) ]=0 (2) 
using the condition SS g /Sg^ = 0. Here, denotes covariant derivative 
and the operator □ is given as 

1 d ( , „„ d 



□ = 



^gdx^V dx h 
with /i, v — 0, 1, 2, 3 and g^ u as metric tensor components. 
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Taking trace of (2) and doing some manipulations, it is obtained that 



+ rR' 



16nGa 



(4) 



R 3(2 + r)(l+r 

with a 7^ to avoid the ghost problem. 

Experimental evidences support spatially homogeneous flat model of 
the universe [12]. So, the line-element, giving geometry of the universe, 
is taken as 

dS 2 = dt 2 -a 2 (t) [dx 2 + dy 2 + dz 2 } (5) 

with a(t) as the scale factor. It gives expansion rate H — a/a. 
In the homogeneous space-time, given by (5), (4) is obtained as 



R + 3-R + 



d-rR 2 1 



16nGa 



+ rR 



(6) 



a R 3(2 + r)(l + r) 
In most of the situations, for example, radiation model, matter- 
dominated model, and accelerated models, we have a(t) as a power- 
law solution yielding R as the power-law function of a(t) . So, it is 
reasonable to take R as 

a n 

with n > being a real number and A being a constant with mass 
dimension 2. 

R, given by (7), satisfies (6), if 

« > /a\ 2 1 r « nr r A „i 

- + (2-rn-n - =— a~ n 8 

a v '\a) 3(2 + r)(l + r)L 8irGna n J v; 

Eq.(8) integrates to 

a\2_ C 1 r_ a nr A _ 

-) ~ _ a 2[3-n(i+r)] ~3( 2 + r )(l + r) I 8nGna{nr + 2 [3 - n(l - r)]) 

^ a-1 (9) 



3n(-n + 2[3-ra(l + r)]) 
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with C being an integration constant. Setting n — 3, (9) is re-written 



as 



a/ 727rGar(2 + r)(l + r) 27(1 + 2r)(2 + r)(l + r)a 3 ' 

(10) 

This is the modified Friedmann equation giving dynamics of the 
universe. The third term on r.h.s (the right hand side) of this equation 
has the form of density for pressureless matter. This term emerges due 
to non-linear term of curvature in the action (1), hence it is termed as 
dark matter density. So, 

SvrG r_A 

3 Pdm 27(l + 2r)(2 + r)(l + r)a 3 ' 1 ] 

The first term, on r.h.s. of (10), emerges spontaneously and second 
term is the combined effect of linear as well as non-linear term of cur- 
vature in the action (1). A very interesting cosmic scenario is obtained 
on using these two terms (first and second) and taking energy density 
p dc as 



Pde B 72nGar(2 + r)(l + r) a3r 



2A 



(12a) 



with 



A lUirGar(2 + r)(l + r)C ^ 



If = Pdc(s) at a = a s such that a 3r = 2A, (12a) looks like 



A 



r 



-a 3r 



a 3r 
l ~2X 



(13) 



Pde-Pdc(s) 7 2^ Gar (2 + r )(l +r )' 

According to WMAP [13], current values of dark matter density and 
dark energy density are Pd m = 0.23p° r and p° e = 0.73p° r with 

p° = (14) 
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H = lOOhkm/MpcSec = 2.32 x l(T 42 /iGeV and h = 0.68. Using these 
observational values in (11) and (13), we obtain 

87TG 



and 



Pdc = Pdc(s) - (Pde(s) - Pdc) ( — ) 



3r 



1 - a 3r /2X 



1 - af/2\ 



(15) 



(16) 



where 



(Pde(s)-p^)(«o r [l-«o72A]) _1 . (17) 



727rGar(2 + r)(l + r) 
Using (11), (15) and (17), Friedmann equation (10) looks like 

1 - a 3r /2X 



(a) = ^ Pdc(s) ~ ^ 



a \ 3r 
a 



+ 0.23# Z 



2 / a 



(18) 



1 - ag r /2A 

If the term, proportional to a -3 , dominates other terms on the r.h.s. 
of (18), it reduces to 



0.23#n 



aoy 
a J 



yielding 

a{t) = a d [l + 0.72H a~ 3/2 (t - t d )f 3 (19) 

which shows decelerated expansion as a < 0. Here a d and t d are con- 
stants. 

When other terms dominate the term proportional to a~ 3 on the 
r.h.s. of (18) and p dc(s ) > p° c , (18) reduces to 

/n\i / n \3r 1 — a 3r /2A 

(j) = <«*>-<«.>(£) - 

which is integrated to 

l + {^I-^-,)}T 1/3 ' Pi> 



1 - ag r /2A 



(20) 



a(t) = a* 
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giving acceleration as a > 0. Here, t* is the time for transition from 
deceleration to acceleration and a* is the corresponding scale factor. 
Moreover, the constant B in (21) is given by 



B 2 



(Pde(s) - Pdc) ( 4"' 1 - «o72A 



(22) 



It is obvious that for t > t*, (18) reduces to (20). So, if Pde(s) < Pd e > 
(20) shows that (^j < for t > t*. It leads to an unphysical situation. 
So, this possibility is rejected and it is concluded that 



Pde(s) > Pdc- 

Further, DE conservation equation is given as 

p dc + 3H(p dc +p dc ) = 0. 

Connecting (16) and (24), it is obtained that 

1 - a 3r /\ 



(23) 



(24) 



Pde = -Pdc + (Pde(s) - Pdc) ( — ) 



3r 



1 - af/2\ 



(25) 



This equation shows that at a = \ l ^ r = a w ,pd c = — Pdc- So, 

3n 



Pdc = "Pdc + (Pdc(s) - Pde) ( — ) 

\a / 



3r 



a/a v 



\ ( do/ dv 



(26) 



This is the EOS for DE, obtained in this model, which is scale factor 
dependent. a s , given above, and a w are related as a s = 2 1 / 3r a w > a 
and a* < a w < ao- 

Now, (16) and (26) yield 



Pde + 3p dc = -2p dc (s) + (Pdc(s) - Pdc) — ' 

vao/ 



Due to inequality (23), (26) and (27) yield 

Pde + Pde > 0, 



- 2(a/a w ) 



3r- 



3r 



(27) 



(28a) 
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Pde + 3p dc < 0, 



(286) 



for a* < a < a w and 



Pde + Pde < 0, (29a) 
Pde + 3p dc < 0, (296) 

for a > a w . 

(28a,b) give quintessence phase of DE, whereas (29a,b) give phan- 
tomce phase of DE. These two phases are divided by 

Pde = -Pde (30) 

at a = a w as given by (27). Thus, we have transition from non-phantom 
to phantom at a = a w > a 

It is given above that at af' = 2A = 2a^ r . So, (18) is obtained as 



(Pde(s)-pD(^) 3r 



2 ( a/a f 



3r 



3r 



^ + 0.23F; 



(31) shows that when a(t) reaches a s , acceleration given by (20) stops 
and deceleration driven by matter resumes. It shows transient acceler- 
ation as obtained in [8]. Sahni [14] had obtained this type of result in 
the context of brane-gravity cosmology earlier. 
(31) also shows that a — at a — a m > a s if 

1 - 2 a m /a s 



(Pde(s) - Pde) (7^) 



3r 



3r 



+ 0.23if 2 (^) =0. (32) 



1 - [a /a s 

It shows that expansion will reach its maximum at a = a m and it will 
begun to contract taking a turn around. 

Thus, a transition from deceleration to acceleration, at some time 
t* in the recent past, is obtained giving a possible explanation for late 
cosmic acceleration [2]. EOS for DE, derived here, depends on the scale 
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factor a(t). Interestingly, it is found that another transition, from 
quintessence phase of DE to phantom phase, takes place at a = a w 
with a* < a w < ao- Characteristics of quintessence DE are different 
from phantom one. Kinetic energy for the latter is negative, whereas 
it is positive for the former. Moreover, Wd c > —1 for the former, 
Wde < — 1 for the latter. So, it indicates that possibly, DE has two 
components (i) quintessence and (ii) phantom. Former dominates when 
a* < a < a w and latter dominates when a w < a < a s . It is found that, 
after acceleration for some time, universe decelerates . The decelerated 
expansion continues till scale factor acquires its maximum at time t m 
and will begin to contract for t > t m . 
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